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Abstract 

In this research the transverse hydrostatic stability of a gravitating Fermi-Dirac plasma under 
a quantizing field is explored. It is revealed that such plasma is magnetically unstable due to the 
Landau electron spin-orbit quantization in all magnetic levels which might explain the gravitational 
compactification in the case of magnetic compact-stars. It is also revealed that the Coulomb inter- 
action plays a fundamental role in the magnetic instability when a critical field is reached for a given 
plasma composition. It is remarked that, the Coulomb interaction, on specific condition depending 
on the field-strength and plasma composition, can suddenly destroy the transverse stability of the 
relativistically degenerate Fermi-Dirac plasma. The consequences of such instability on mass-limit 
of magnetic compact-stars and the neutron star crust with possible relation to gamma-ray bursts 
is numerically examined. 

PACS numbers: 
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I. INTRODUCTION 



Astrophysical compact plasma is a unique environment in which newtonian gravity, spe- 
cial relativity, quantum mechanics and electromagnetism all play equal parts. Perhaps the 
most distinguished exhibition of large-scale classical gravity and quantum mechanics col- 
laboration is the effect of relativistic degeneracy of Fermi-Dirac electrons leading to the 
famous Chandrasekhar mass-limit on compact-stars called white-dwarfs It has been 
shown that the equation of state Q as well as nonlinear wave dynamics jsj] of a degenerate 
plasma is affected fundamentally by quantum degeneracy when the de Broglie thermal wave- 
ength Xb = /i/(27rm e /csT) 1 / 2 becomes comparable or higher than inter-fermion distances 
4J. There are other effects such as Thomas- Fermi crystallization, Coulomb interaction, 
electron exchange and ion correlation effects which might contribute to stability of quantum 
plasmas. However, it has been shown [5[ that these effects give-rise to correction-terms of 
order the fine-structure constant or even smaller to the hydrostatic stability of gravitating 
relativistically degenerate plasmas slightly altering the mass-radius relation only for high 
atomic-number composed compact stars 6(. 

One of the greatest challenges in astrophysical sciences is to understand puzzling presence 



and possible role of magnetism in stellar birth 
physical processes such as gamma-ray bursts 
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collapse, js|, evolution 10, [11] and 
13|. There have been reports of white- 

15fl • Extensive 



dwarfs and neutron stars with associated ultra-high magnetic fields [ijj 
past studies 16K18I] reveal that the magnetism can lead to quantization of thermodynamic 
quantities and electronic density of state in Fermi-Dirac gas. The magnetic instability has 
been shown to dominate due to anisotropic stress in ground-state Landau-level 19j. The 
oscillatory feature of spin-orbit magnetization in super-dense plasma may also lead to a 
metastable self-sustained ferromagnetism (called Landau orbital ferromagnetism or LOFER) 



with field strength as high as 10 8 G and 10 13 G for white-dwarf and neutron-star density ranges 
20, [2lJ. However, the LOFER state has been shown to disappear at temperatures estimated 
for such compact-stars 22|. Recently, it has been shown that the magnetic transverse 
collapse due to pressure anisotropy in magnetized neutron-star 23( and LOFER plasma 24 1 
is possible. 

One of important astrophysical motivation towards the study of matter under extreme 
conditions such as in strong magnetic-field and relativistic degeneracy arises due to im- 
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portance of understanding the neutron star crust which play fundamental role in many 

n 

astrophysical processes and observed phenomena such as gamma-ray bursts [25]. The study 
of matter under stron g m agnetic field has also implication for the strong laser-matter in- 



teraction experiments 
for gamma-ray bursts 



Many models has been proposed for the underlying mechanism 
27fl . However, the inner engine for such gigantic energy release is far 



from clear. More recently, a model of resonant excitation in the neutron star (NS) crust has 



been proposed [28j . It has been shown that the crust-core interaction may lead to cracks in 
the crust leading to energy release of order 10 46 ergs. However, the crust tension in this and 
many other previous models rely on the external mergers such as that of NS-NS or NS-black 
hole. In the present paper we show that the instability and shatter of a neutron-star crust is 
possible due to Landau quantization accompanied by Coulomb interactions without a need 
for external cause. 



II. STABILITY CRITERION 



Considering a uniformly magnetized homogenous gravitating Fermi-Dirac quasi-neutral 

plasma, the hydrostatic stability is defined through the force equilibrium, F in = F out on an 

spherical shell element, dr, may be written as 

1 dP _ GM(r) 
p dr r 2 ' 



(1) 



which at the quasi-neutral state, i.e. n e ~ rij, the local center of mass plasma density will be 
p = rriini + m e n e ~ m,in e . The local pressure of relativistically-degenerate electrons compose 
almost the total pressure which can be anisotropic, P\\ = P e = P± + TB, in the presence of a 
uniform magnetic field, in which T is the magnetization. The electron pressure and density 
are Landau-quantized and at the zero-temperature approximation take the form 16j 



p(e Fe , 7) = p c 7 J2 ( 2 - ko) V e Fe -2/7-I, 

1=0 

P\\(e Fe ,j) = \p c ^m e c 2 ( 2 - Si fi ) \ e FeA /e 

1=0 
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Fe 



2/7 



-(I + 2/7) In 
r(e Fe ,7) = 5- 1 



tFe+y/e 2 Fe -2h-l 



<3P| |(eF e ,7) 

d-y 



(2) 



where, p c = m i vr? e c' /2ir 2 h 3 , 7 = B/B c with B c = m 2 e c 3 /eh ~ 4.414 x 10 13 G. The parameter 



ep e — vl + x 2 is the normalized relativistic Fermi-energy which is related to the relativity 



parameter, x = pF e jm e c, the normalized Fermi electron-momentum. The maximum oc- 
cupied Landau-level, Z m , is defined as l m = (e 2 Fe — l)/2 r y > Z and the parameter 5i t o takes 
values of {1,0}, respectively, for the quantum-limit (/ = 0) and for other /-values. The 
corresponding Chandrasekhar equation-of-state may be obtained from Eq. (T2]) as / — > oo. 
The transverse hydrostatic stability can be written as 



d ( r 2 d[P {l (x,j)-BT(x^)] 



4nr 2 G dr \ p(x, 7) dr 



+ p(x, 7 ) = 0, (3) 



On the other hand, the equation of state of a magnetized Fermi gas can be evaluated 
analytically in terms of Hurwitz functions [29] in the following form 

p(*,7)=p c (27) 3/2 #-i/ 2 (g), 

P || (x, 7 ) = ^p c (27) 5 / 2 #^¥i|^ 

H p (q) = h{p, {q}) - h(p, q + 1)- \q~ p , 

00 

hip,q) = E {n + q)~ p . 

where hip, {q}) is the Hurwitz zeta-function of order p with the fractional part of q 
as argument and c s = c\J m e /rrii is the wave-speed. Using the definitions p(e j p e ,7) = 
c~ 2 dP\\ iep e , ^)/deFe and r(ei? e ,7) — ^ 1 ^|( e -Fe,7)/^7 for a homogenous magnetized 
plasma case, we have 



1 <9[P||0r,7)-r(x,7)£?] _ c 2 s xir) 



P{x,l) dr Jl + X ( r ) 2 



' _ 7 dpjx,^) 
p(x,7) ^7 



which, after some algebra, leads to the general stability relation below for the uniformly 
magnetized gravitating Fermi-Dirac plasma 



c 2 d I r 2 x(r) 



Airr 2 G dr 



\Jl + x(r)' 



_ 7 dpjx,^) 
p(x,7) d-f 



dxir) 
dr 



+ p(x, 7 ) = 0. (6) 



The parameter £7(2,7) = 7<91np(x, "f)/d^ defines the stable solution to Eq. (jSJ). Figure 1 
shows the variation of mass-density (along with the monotonic and oscillatory components) 
and stability parameter n(ar, 7) in terms of the relativity parameter, x, indicating that 
the stability solution can become unstable for some discrete x-ranges for a given magnetic 
field parameter, 7. The linear part of Fig. 1(a) corresponds to the quantum-limit (Z = 0) 
for which p(x, r {) = jx. On the other hand, at the weak quantization limit (Z = 00) the 



magnetization vanishes and we have p = x 3 , where, the plasma becomes magnetically stable 
and isotropic in all directions. The corresponding Landau-level is shown in Fig. 1(b) to 
be critically-unstable (0(x, 7) = 1). Landau-levels along with the corresponding quantized 
stability regions are depicted in Fig. 2. The white regions in Fig. 2(b) indicate the stability 
voids for which tt(x, 7) < 1. For the quantum-limit, Q(x, 7) = 1, (the largest white-region 
to the left in Fig. 2(b)) the plasma can free-fall under its own gravity, while, for other 
Landau-levels the total pressure becomes negative and the plasma collapse may occur even 
in the absence of gravity. 



III. COULOMB EFFECT 



It has been shown that the effects such Thomas-Fermi screening, Coulomb attraction, 
exchange and correlation effects can affect the stability of the relativistically degenerate 
plasma jsj. The former three effects are of the order of fine-structure constant, while the 
correlation effect is much smaller. Here, we consider the Coulomb effect which can be 
dominant particularly in the presence of a magnetic field, since, the degeneracy of the plasma 



can be dominated by the new bound-states introduced by Landau quantization [25j]. The 
attractive Coulomb interaction introduces another negative term to the pressure of plasma, 
as below p] 



3tt 



1/3 



(7) 



„ . > 187r 2 m 4 e c 5 fa 5 z 2 l 3 \ 
Pc(Xn) = ~ U^J' ^ a [S 1P (x^\ 

where, the parameters, a = e 2 /hc and z are the fine-structure constant and the atomic 
number, respectively. Hence, in the presence of Coulomb effect the stability parameter can 
be written as 

A _ 7 <9p(x,7) | g>/j + x 2 z 2/3 7 4/3 dp(x,7) _ 647r 2 / 3 q 
lX,7 '" j "p(*,7) d 1 + xpix^f* dx ' "-45X3VS- ^ 



where, we have used the identity 

1 dPc dp(x, 7) dx dtFe 



7 c 



(9) 



p(x, 7) dp(x, 7) dx dtFe dr 
The stability regions for two different plasma compositions, namely carbon and iron, are 
shown in Fig. 3. It is remarked that for a given critical- value of fractional field-parameter, 
7, the stability is shattered. Surprisingly, the Coulomb stability-shatter for iron composition 



is observed to coincide with critical-field value, B ~ B c . Evidently, this effect can lead to 
shattering in the external iron-rich crust of a critically-magnetized neutron star and allow 
to a enormous energy release such as gamma-ray bursts. In a recent paper Tsang et al. [| 
have shown that such energy release is possible due to a resonant shatter in neutron-star 
crust. 



IV. FINITE-TEMPERATURE EFFECT 

The finite-temperature effects might be important for some white-dwarfs with typical 
temperature of order 10 6 Kelvins. The thermodynamic quantities of finite-temperature may 



be calculated from the zero-temperature counterparts, through the relation 29| 



Q(0,fi + ^r) k B T 
4cosh 2 (£/2) T m e c 2 



Q( T ^)= I ,...^2^^ d Z, r= ZT^' ( 10 ) 



where, \x is the chemical potential and Q can be an arbitrary thermodynamic quan- 
tity. For instance the density p(/i, r, 7) can be obtained from p(x', / y) with argument 
x = \J x 2 + 2/ixr + /i 2 r 2 , which ignoring the p 2 r 2 for small values of r leads to 

p(T '"' 7) J i^W) 1/2 I - »r~J * ( ' 

—00 

For a comparison Fig. 4 shows the effect of finite-temperature (ignoring the Coulomb interac- 
tion effect) on plasma-density. It is observed that the oscillatory part of the thermodynamic 
quantities introduced by h(p, {q}) in Eq. is smoothed in finite-temperature case, how- 
ever, the magnitude is not affected significantly. Due to such effect one expects that the 
stability quantization in Fig. 2(b) to be smeared out at high temperatures. To evaluate 
the finite-temperature effect on plasma instability one may consider only the monotonic 
part of density function in Eq. (j^). Figure 5 shows the stability regions for the above- 
mentioned compositions taking into account only the monotonic part of density function. It 
is remarked that, the destabilization effect persists even with monotonic variation of ther- 
modynamic quantities and is not specific to the quantized case. This is an indication of the 
fact that the Coulomb instability can even be present at finite-temperature case. Also, the 
values of field-parameters for which the instability starts in Fig. 5 are of the same order of 
magnitude in Fig. 3. 



V. THE MASS-LIMIT 



To consider the aforementioned instability on mass-limit we will consider only the mono 



tonic variation and calculate based on the method given in Ref. 



30] . The single differential 



equation involving the stellar profile transverse to the magnetic-field is of the form 

^ f x& 11 " M ^} + (27)3/2ff - i/2 M " °' (12) 



where, the stability parameter Q(x, 7, z) is defined in Eq. and n h = ^Jhc/G /rrii ~ 1.3 x 
10 19 is the dimensionless hierarchy-number defined based on three fundamental quantum, 
relativity and classical constants, namely, the scaled Plank-constant, h, the speed of light 
in vacuum, c and the gravitational constant, G. It is clearly evident that, the hydrostatic 
stability of the configuration given by Eq. ( Tl2l is lost when the term in the brackets in this 
equation becomes negative or vanishes. We now move-on to evaluate the stable configuration 
for given values of relativity and field parameters. 

Each stable configuration is found by integrating Eq. ([12]) from the center (r = 0) 
outwards until x(r) vanishes for some r = r s which would be the surface (of the star). The 
initial condition, however, is set to x(r = 0) = x c and dxjdrir = 0) = such that for 
every given value of the relativity parameter at the center (x c ) there will be a distinct stellar 
configuration. It is convenient to change the Eq. (fl2l) to the normalized form by change of 
variables y = x/x c and £ = r/r c , where, r c = \/37inhX c /4:X c , so that at the center (£ = 0) we 
will have y = 1. The Eq. f|T2|) . then in dimensionless form, reads as 

1 d j x c ey(Q n m „ dy(Q 1 (2 7 ) 3 / 2 ( xM£f \ n n ,. 



l + x^(0 2 ' ' * J 

It is remarked (making use of Eq. (jSj)) that, in the field- free limit, the above relation reduces 
to the famous Lane-Emden equation of index 3 leading to the well-known Chandrasekhar 
mass-limit of Mch — 1.43217Mg (Ms being the mass of the sun). A standard integration 
algorithm may be used to find, for instance, the value of £ = £ s at which ?/(£) vanishes (the 
surface of star) from which the radius and the mass of the stable configuration (star) can 
be calculated via the following relations 

r s = r c £ s , M = 4vr £ p(r)r 2 dr = ^m^f H_ 1/2 (^f-) ^ (14) 



In the next section, for numerical scheme we will use the unique scalings 7ihX c ~ 0.787Re 
(Re being the radius of Earth) and min\ ~ 1. 849 Ms (Ms being the mass of Sun), which are 
conventionally employed in literature. Figure 6 shows the solution of y(£) of Eq. ( |14p for 
different values of fractional field-parameter, 7, for dense-core configuration (x c ~ 00) for 
uniform iron composition. It is evident that the mass-limit varies according to the change 
in magnetic field strength. However, the configuration becomes unstable (dashed curves) 
for values of 7 for which the stability is shattered (e.g. see Fig. 3). The composition 
dependent of the normalized mass-limit id depicted for three different compositions in Fig. 
7 making a bold conclusion about the well-known Chandrasekhar mass-limit. It should be 
noted that although the field strength considered here may be vary high compared to those 
measured for white-dwarfs surface, these values can be realized for neutron stars and even 
for white-dwarf cores 31]. 
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FIGURE CAPTIONS 



Figure- 1 

Variation of the normalized plasma mass-density (along with oscillatory and monotonic 
parts) and stability parameter, Q(x, 7), with respect to the fractional Fermi-momentum, x 
(the relativity parameter) for a given magnetic field parameter, 7 = 0.3. 

Figure-2 

Regions of the Landau magnetic-levels and the corresponding regions of stability voids 
(unstable regions in white) for a uniformly magnetized zero-temperature Fermi-Dirac plasma 
caused by the spin-orbit quantization. 

Figure-3 

Coulomb stability-shatter for two different plasma composition, namely, the carbon and 
iron. The Coulomb stability-shatter starts at a given critical magnetic-field value (7 ~ 2.5 
for carbon and 7 ~ 1 for iron plasma composition) 

Figure-4 

Effect of finite-temperature (T ~ 5.9 x 10 6 fT) on Fermi-Dirac plasma mass-density quan- 
tization for a chemical potential value, \i = 100. The oscillatory nature of the plasma 
quantities are smeared-out as the plasma temperature is increased. 

Figure-5 

Effect of the Coulomb instability in the absence of oscillatory part of plasma mass-density. 
The stability region (colored regions) up to very large values of relativity parameter, x, 
is bounded by critical field-parameter depending on the plasma composition. The field 
parameter values in these plots are slightly higher but consistent with the ones in Fig. 
3. Also the left stability branches in these plots are due to approximation (ignoring the 
oscillatory part of the mass-density) which are not present in Fig. 3. 

Figure-6 
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The solution to for dense-core iron composition plasma with different values of 
magnetic field strengths. The solid and dashed curves indicate the stable and unstable 
configuration, respectively. The thickness of curves increases according to increase in the 
value of 7. 

Figure-7 

Effect of plasma composition and the strength of magnetic field on the compact-star mass- 
limit. The cases of z — 2, 12, 56 shown respectively with curves of rectangular, filled- circles 
and halo- circles symbols. 
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